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Abstract 



In this paper we propose an interpretation of the time series of 
the RMS overlap differences of the out-of-plane part N of the orbit of 
the Martian polar orbiter Mars Global Surveyor (MGS) over a time 
span AP of about 5 years (14 November 1999-14 January 2005) in 
terms of the general relativistic gravitomagnetic Lense-Thirring effect. 
It turns out that the average of such a time series over AP, normalized 
to the predicted Lense-Thirring out-of-plane mean shift over the same 
time span, is n = 1.0018 ± 0.0053. A linear fit of the considered time 
series over AP yields a normalized slope 1.1 ± 0.3 strengthening our 
interpretation of the determined MGS out-of-plane features of motion 
as due to the gravitomagnetic force of Mars. 
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1 The gravitomagnetic Lense-Thirring effect 

In the framework of the weak-field and slow-motion linearized approximation 
of general relativity [31], the gravitomagentic Lense-Thirring effect consists 
of tiny secular precessions of the longitude of the ascending node O and 
the argument of pericentre uj of the orbit of a test particle freely moving 
in the gravitational field of a slowly rotating body of mass M and angular 
momentum J [27, 2, 9, 46, 1, 13] 
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where G is the Newtonian gravitational constant, c is the speed of light in 
vacuum, a, e and i are the semimajor axis, the eccentricity and the incli- 
nation of the orbital plane to the central body's equator, respectively. It 
maybe of interest to recall that de Sitter in [10] worked out, for the first 
time, the gravitomagnetic pericentre precession in the case of equatorial or- 
bits (i = deg). In general, the gravitomagnetic force affects the in-plane 
transverse and the out-of-plane cross-track portions of the orbit of a test 
particle. Indeed, the shifts along the radial R, transverse T and normal 1 
N directions can be, in general, expressed in terms of the Keplerian orbital 
elements as [3] 
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where A4 is the mean anomaly The Lense-Thirring shifts are, thus 
( AR LT = 0, 
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For an equatorial orbit there is only the transverse shift ATlt, while for a 
polar (i = 90 deg) satellite only the out-of-plane shift ATVlt is present. This 
fact can be more intuitively grasped in the following way. In general, a time- 
variation of the argument of pericentre to, which is an angle entirely residing 
in the orbital plane delimited by the line of the nodes and the apsidal line, 
can be obtained by subtracting the projection of the nodal motion to the 
orbital plane, i.e. cositi, from the precession of the longitude of the peri- 
centre w. For a polar orbital geometry the projection of the node's motion 
to the orbital plane is zero. According to the well-known Gauss variation 
equations [43], the remaining shift induced by a small perturbing accelera- 
tion A is due to its in-plane radial Ar and transverse At components. The 



1 It is customarily denned to be positive along the orbital angular momentum. 
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Lense-Thirring effect is due to a Lorentz-like acceleration [28] 
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the gravitomagnetic field [28] 

_ G[3r(r-J)-r 2 J] 

for i = 90 deg lies in the test particle's orbital plane which contains also J, 
so that the net gravitomagnetic force is entirely normal to the orbital plane. 

Up to now, many attempts to test gravitomagnetism were implemented 
in the weak-field and slow-motion scenarios of the gravitational fields of some 
Solar System's bodies. Lense and Thirring considered the natural satellites 
of Jupiter and Saturn [27], but the possibility of using the jovian galilean 
moons is not yet viable [22] . In regard to the Sun and its rocky planets [20] , 
the predictions of general relativity for the Lense-Thirring perihelia preces- 
sions were, in fact, found in agreement with the estimated corrections to the 
secular rates of the perihelia of the inner planets [38] , but the errors are still 
large. As far as the Earth is concerned, a test of the Lense-Thirring effect 
on the orbit of a test particle was reported in [6]. Their authors analyzed 
the data of the laser-ranged LAGEOS and LAGEOS II satellites by using an 
observable independently proposed in [41, 42] and [21] involving the nodes of 
the LAGEOS satellites. The idea of using the LAGEOS satellite and, more 
generally, the Satellite Laser Ranging (SLR) technique to measure the Lense- 
Thirring effect in the terrestrial gravitational field with the existing artificial 
satellites was put forth for the first time by Cugusi and Proverbio in 1977- 
1978 [8, 9]; tests actually started in 1996 with LAGEOS and LAGEOS II [5]. 
Although the outcome of the test in [6] is in agreement with the predictions 
of general relativity, the realistic evaluation of the obtained total accuracy 
raised a debate [15, 16, 17, 19, 7, 29]. The GP-B spacecraft [11, 12] was 
launched in April 2004 with the aim of measuring another gravitomagnetic 
effect, i.e. the Schiff precession [39, 44] of the spins of four superconduct- 
ing gyroscopes carried onboard. Preliminary results, released in April 2007, 
yield a measurement uncertainty of the order of 256 — 128%; the final results 
should be announced within the end of 2007, with a goal of an about 13% 
accuracy, (http://einstein.stanford.edu/). Recently, a preliminary 6% LT 
test on the orbit of the Mars Global Surveyor (MGS) spacecraft (see Table 1 
for its orbital parameters) in the gravitational field of Mars was reported in 
[18]; indeed, the predictions of general relativity are able to accommodate, 
on average, about 94% of the Root-Mean-Square (RMS) overlap differences 
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Table 1: Orbital parameters of MGS [24]; a is the semimajor axis, i is the 
inclination to the Mars'equator and e is the eccentricity. 



a (km) i (deg) e 
3796 92.86 0.0085 



in the out-of-plane part of the MGS orbit over 5 years. Finally, it must 
be mentioned that, according to Nordtvedt [35, 37], we would already have 
evidences of the existence of the gravitomagnetic interaction. Indeed, the 
multi-decade analysis of the Moon'orbit by means of the Lunar Laser Rang- 
ing (LLR) technique yields a comprehensive test of the various parts of order 
0(c~ 2 ) of the post-Newtonian equation of motion. The existence of gravit- 
omagnetism as predicted by general relativity would, thus, be inferred with 
high precision from the reconstruction of certain periodic terms in the lunar 
motion. A 0.1% accuracy is claimed in [33], but, according to Kopeikin [25], 
LLR would not be capable to detect gravitomagnetic effects; Murphy et al. 
replied in [34]. Other comments on [33] are in [4]. Also the radial motion of 
the LAGEOS satellite would yield another confirmation of the existence of 
the gravitomagnetic interaction [36]. 

In this paper we perform a test of the Lense-Thirring effect by improving 
the accuracy of the previous one with MGS [18] by means of a different data 
set, a refinement of the overall calculation and a realistic and conservative 
assessment of the error estimate. The details of the adopted strategy can 
be found in [18, 24]; here we will briefly recall them. 

2 Using the MGS data in the gravitational field of 
Mars 

In [24] six years of MGS Doppler and range tracking data and three years 
of Mars Odyssey Doppler and range tracking data were analyzed in order 
to obtain information about several features of the Mars gravity field sum- 
marized in the global solution MGS95J. As a by-product, also the orbit of 
MGS was determined with great accuracy In Figure 3 of [24] the time series 
of the RMS overlap differences of the radial, transverse and normal portions 
of its orbit over the entire mission (4 February 1999—14 January 2005) are 
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shown; their averages over such a time span amount to [24] 

' (AR diS ) = 15 cm, 

< (ATdiff) = 1.5 m, (6) 

, (AJVdiff) = 1.6 m. 

The RMS overlap differences are customarily used in satellite geodesy as 
indicators of the overall orbit accuracy [48], concisely accounting not only 
for the measurement errors, but also for all the systematic bias due to the 
mismodelled/unmodelled forces acting on a spacecraft [30]. Errors common 
to adjacent arcs are, by construction, cancelled; it is not so for effects grow- 
ing in time like the one we are interested in here. Thus, any analytical 
calculation of possible systematic errors must ultimately cope with that is 
actually shown by the RMS overlap differences. It is important to note that 
in the suite of dynamical force models used to reduce the gathered observa- 
tions general relativity is, in fact, present, but only with the gravitoelectric 
Schwarzschild-like part: the gravitomagnetic one is not included, so that, in 
general, it is fully accounted for by the RMS orbit overlap differences. From 
eq. (3) it turns out that only the out-of-plane portion of the orbit of a polar 
satellite like MGS is affected by gravitomagnetism via the Lense-Thirring 
effect: from eq. (1) and eq. (3) the mean gravitomagnetic normal shift, 
averaged over a temporal interval AP, is 2 

Moreover, no empirical out-of-plane accelerations were fitted in the data re- 
duction process; this was only done for the along-track component of MGS 
orbit to account for other non-gravitational forces like thermal, solar radi- 
ation pressure and drag mismodelling. Thus, it is unlikely that the out-of- 
plane Lense-Thirring signal was removed from the data in the least-square 
procedure. Moreover, if the relativistic signature was removed so that the 
determined out-of-plane RMS overlap differences were only (or mainly) due 
to other causes like mismodelling or unmodelling in the non-gravitational 
forces, it is difficult to understand why the along-track RMS overlap differ- 
ences (middle panel of Figure 3 of [24]) have almost the same magnitude, 

2 A 1/2 factor appears as a result of the time-average of a linear signal like the gravit- 
omagnetic one over a time interval AP. 
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Table 2: Relevant parameters of Mars used in the text. The parameter 
a = C/MTZ 2 is the normalized polar moment of inertia, 1Z is the reference 
Mars radius in MGS95J model, while V is the Martian sidereal period of 
rotation. From such values and from J = aMTZ 2 2ir/V it can be obtained 
J = (1.92 ± 0.01) x 10 32 kg m 2 s _1 for the Martian angular momentum. 

GM (10 a m 3 s~ 2 ) Tl (km) a V (d) 

42828.37440 ± 0.00028 [24] 3396.0 [24] 0.3654 ± 0.0008 [24] 1.02595675 [45] 



since the along-track component of the MGS orbit is much more affected 
by the non-gravitational accelerations (e.g. the atmospheric drag) than the 
out-of-plane one. Finally, the following consideration should further enforce 
our conclusion. The improved modelling of the non-gravitational forces act- 
ing on MGS introduced in [24] with respect to previous works [26, 51] in 
which the Lense-Thirring effect was not modelled as well, only affected in 
a relevant way the along-track RMS overlap differences (a factor 10 better 
than [26, 51]), not the normal ones (just a factor 2 better than [26, 51]). 

By using a time span AP of slightly more than 5 years (14 Novem- 
ber 1999—14 January 2005), so to remove the first months of the mission 
more affected by orbital maneuvers, high gain antenna deployment and non- 
gravitational perturbations, it is possible to obtain 

(AiV diff ) = 1.613 m. (8) 

Let us, now, compare the result of eq. (8) with the mean gravitomagnetic 
shift predicted by eq. (7) over the same time span; according to our hypoth- 
esis, the Lense-Thirring force of Mars is responsible for the averaged normal 
signature of eq. (8). The figures of Table 2 for the Mars'relevant parameters 
and those of Table 3 for the fundamental constants entering the calculation 
yield 

(AiV diff ) 



/^mcas — 



= 1.0018 ±0.0053; (9) 

AP 

according to our hypothesis, fJ, e xpected = 1- The error was conservatively 
evaluated as 



(AiV LT ) 
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Table 3: Values used 


for the denning, primary 


and derived 


constants 


(http: / / ssd.jpl.nasa.gov 


/?constants#ref). 






constant 


numerical value 


units 


reference 


c 


299792458 


m s _1 


[32] 


G 


(6.6742 ± 0.0010) x 10~ n 


kg -1 m 3 s~ 2 


[32] 


1 mean sidereal day 


86164.09054 


s 


[47] 


1 mean sidereal year 


365.25636 


d 


[47] 



with 



f 8n\j < fi (^) = 0.0052, 
8h\g < H (7?) = 0.0001, 



(11) 



I MU<2 /U (^)=O(10- 8 ). 

As can be noted, the major source of bias is represented by the angular 
momentum of Mars. Its uncertainty was evaluated as 



SJ <SJ\ n + SJ\ a + 6J\ M , 



(12) 



with 



' sj\ n <2j(^) =7xl0 29 k 



g m 2 s 1 , 



SJ\ a <J(^)= 4 x 10 29 kg m 2 s" 



(13) 



SJm < J 



8{GM) 
GM 



0(1O 24 ) kg m 2 s" 



We linearly added all the errors in a conservative way; in regard to a and 
GM this is particularly important because they are fit-for parameters of 
the global solution and, consequently, are certainly affected by correlations 
which were not explicitly released in [24]. It must also be noted that the 
uncertainties quoted in Table 2 are realistic, not the mere formal, statistical 
errors; indeed, the quoted uncertainty in a is, e.g., 15 times the formal one 
[24]. Also the formal error in GM was re-scaled [24]. We also considered 
an uncertainty in 7Z by assuming as a realistic measure of it the difference 
51Z = 6080 m between the reference MGS95J value for the equatorial radius 
quoted in Table 2 from [24] and the value 3389.92 km released in [50], instead 
of the error of 0.04 km quoted in [50]. In regard to the error in a, we assumed 
from eq. (6) 5a ~ 15 cm in view of the very small eccentricity of the MGS 
orbit. 
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A further confirmation of the validity of our interpretation of the MGS 
out-of-plane motion in terms of the gravitomagnetic force of Mars comes 
from a linear 3 fit of the time series of the RMS overlap differences over AP. 
Indeed, the estimated slope is —0.69 ± 0.20 m yr -1 , while general relativity 
predicts a Lense-Thirring out-of-plane rate of 0.62 m yr _1 , with the normal 
direction defined positive along the orbital angular momentum. The minus 
sign of the fit is due to the fact that in the analysis of [24] the MGS normal 
direction is, instead, defined positive in the opposite direction (Konopliv A, 
private communication 2007). Incidentally, let us note that should we as- 
sume such a fit as indicator of the existence of the Lense-Thirring effect, the 
predictions of general relativity would be confirmed within the experimental 
error which, however, would be larger than the one obtained by taking the 
mean of the time series of N over AP. 

Another important remark concerning our measurement is the follow- 
ing one: the interpretation of the MGS orbital features presented in this 
paper and their production from the collected data were made by distinct 
and independent subjects in different times; moreover, the authors of [24] 
had not the Lense-Thirring effect in mind at all when they performed their 
analysis. Thus, the possibility of any sort of (un)conscious gimmicks, tricks 
or tweaking the data because the expected answer was known in advance 
[49] is a priori excluded in this test. 

3 Conclusions 

We conclude by summarizing our main result: the ratio \x of the mean of 
the time series of the cross-track N RMS overlap differences of MGS from 
14 November 1999 to 14 January 2005 to the mean of the Lense-Thirring 
shift over the same time interval amounts to [i = 1.0018 ± 0.0053. Our 
interpretation of the features of the out-of-plane motion of MGS in terms of 
the gravitomagnetic force of Mars is enforced by a linear fit of the N time 
series, although the error would be, in this case, larger. 

It must be noted that the present analysis is based on the so far pro- 
cessed data in [24]; other approaches which could fruitfully be used in the 
future may consist in re-processing the entire MGS data set including in 
the dynamical force models the gravitomagnetic force as well and checking 
if the average of the N time series gets smaller, and/or including a further 
solve-for parameter in the global least-squares solution accounting for the 

3 It may be just the case to note that, since the plots in Fig. 3 of [24] are semi- 
logarithmic, one should not visually look for a straight line in them. 
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Lense-Thirring effect itself. 

Finally, let us note that, although we lost contacts with MGS since early 
November 2006, more independent tests could be, in principle, performed 
in the near future when a sufficient amount of data from the Odyssey and, 
especially, Mars Reconnaissance Orbiter (MRO) orbiters will be available 
and processed. 
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